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Abstract We describe an extension of the Feyn Rules 
package dedicated to the automatic generation of the 
mass spectrum associated with any Lagrangian-based 
quantum held theory. After introducing a simplified 
way to implement particle mixings, we present a new 
class of FeynRules functions allowing both for the 
analytical computation of all the model mass matri- 
ces and for the generation of a CH — h package, dubbed 
ASperGe. This program can then be further employed 
for a numerical evaluation of the rotation matrices nec- 
essary to diagonalize the field basis. We illustrate these 
features in the context of the Two-Higgs-Doublet Model, 
the Minimal Left-Right Symmetric Standard Model and 
the Minimal Supersymmetric Standard Model. 



1 Introduction 

Although the Standard Model of particle physics is very 
well verified empirically at the current accessible en- 
ergies, numerous extensions to its Lagrangian are pro- 
posed. These extensions describe new or alternative fun- 
damental interactions that typically accommodate pos- 
sible new physics phenomena at higher energies as well 
as at the current collider energies. In this top-down ap- 
proach, the phenomenology of the proposed extensions 
is to be confronted with experimental observations. In 
order to obtain the mass spectrum of any new physics 
model reflected by its Lagrangian, the mixing matri- 
ces of the gauge eigenstates into the mass eigenstates 
are needed. An automated mass spectrum generator, 
ASperGe 1 is developed within the framework of the 
FeynRules program [1,2,3,4,5] to determine the mix- 



ing matrices numerically. This allows a study of the di- 
rect relation between the parameters of any new physics 
model and the observable masses of the fundamental 
particles. 

This paper describes in Section 2 and 3 the relevant 
parts of the FeynRules program to introduce the new 
ASperGe package, detailed in Section 4. To illustrate 
its application several examples are presented in Sec- 
tion 5. 



2 The FeynRules package 

The program FeynRules [1,2,3,4,5] is a Mathemat- 
ICA 2 package that allows for the automated extrac- 
tion of Feynman rules from any Lagrangian describing 
the dynamics of a perturbative quantum field theory. 
The Feynman rules, together with general information 
such as the definitions of the model particles or of the 
Lagrangian parameters, can subsequently be exported 
by means of several translation interfaces to matrix- 
element generators. Up to now, interfaces to CompHep 
and CalcHep [6,7,8,9], FeynArts and FormCalc 
[10,11,12,13], MadGraph and MadEvent [14,15,16, 
17,18], Sherpa [19,20] and Whizard [21,22] have been 
developed. 

In addition, any model can also be converted to a 
Python library containing classes and objects repre- 
senting particles, parameters and vertices. This format 
is dubbed the Universal FeynRules Output (UFO) 
format [23] and is appropriate to address the imple- 
mentation of any high-energy physics model into com- 
putational tools. Its strength lies in its agnosticism with 
respect to the allowed Lorentz and/or color structures 



x The acronym ASperGe stands for Automated Spectrum 
Generation. 



2 Mathematica is a registered trademark of Wolfram Research, 
Inc. 
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appearing in the Lagrangian, in contrast to any other 
more conventional model format for which restrictions 
are imposed. Presently, the UFO is used by Aloha [24] . 
Mad Analysis 5 [25], and MadGraph 5 [18], and will 
be used, in the future, by GoSam [26,27] and Her- 
wig++ [28]. 

The FeynRules model structure extends the for- 
mat employed in FeynArts [11] so that particles, pa- 
rameters and gauge groups are now all defined in a 
similar fashion. Following the FeynArts conventions, 
particles are collected into classes describing multiplets 
having exactly the same quantum numbers, but possi- 
bly different masses. Each of these classes is defined as a 
set of Mathematica replacement rules referring to its 
properties. For example, the three vector fields W, as- 
sociated to the SU (2) l gauge subgroup of the Standard 
Model could be declared as 

V[l] == { 

ClassName -> Wi, 

Unphysical -> True, 

Self Conjugate -> True, 

Indices -> { Index [SU2W] }, 

Flavorlndex -> SU2W, 

Definitions -> { 
Wi [mu_ , 1] -> (Wbar [mu] +W [mu] ) /Sqrt [2] , 
Wi [mu. , 2] -> (Wbar [mu] -W [mu] ) / (I*Sqrt [2] ) , 
Wi [mu_,3] ->cw Z [mu] + sw A [mu] } 

> 

This set of Mathematica rules defines a vector field 
(V[l]) represented by the symbol Wi (its ClassName) 
and carrying a flavor index SU2W associated with the 
adjoint gauge index of SU(2)l. As this field is declared 
as unphysical 3 (Unphysical->True), it must be linked 
to one or several of the mass eigenstates of the model 
by means of appropriate mixing relations. As illustrated 
in the example, these relations are passed through the 
attribute Definitions of the particle class. They can 
either be purely numerical, as for the W\ and W2 bosons 
that are rotated to the charged W + and W~ bosons, 
or depend on some of the model parameters, as for the 
W3 field which is re-expressed in terms of the photon A 
and the Z boson through a relation depending on the 
sine and cosine of the electroweak mixing angle (sw and 
cw). 

In contrast, declaring physical particles requires the 
implementation of extra information such as their masses 
(Mass), widths (Width) and Particle Data Group codes 
(PDG) [29]. The Z-boson field introduced above could 
hence be declared as 



V[2] 



{ 



ClassName 

Self Conjugate 

Mass 

Width 

PDG 



-> Z, 
-> True, 

-> {MZ, 91.1876}, 
-> {WZ, 2.4952}, 
-> 23 



} 



The declaration of the model parameters and gauge 
group is similar and based on dedicated classes with 
their own set of attributes. Since only the particle class 
properties introduced above are sufficient for the under- 
standing of the present work, we omit any further detail 
and refer the reader to Refs. [1,4] for more information 
on particle, parameter and gauge group implementation 
in FeynRules. 

The last key ingredient to achieve a model imple- 
mentation consists of its Lagrangian. It is provided us- 
ing standard Mathematica commands, augmented by 
some special symbols representing objects such as Dirac 
matrices, vector field strength tensors or covariant de- 
rivatives. The user has then the possibility to perform 
basic checks on the Lagrangian, such as verifying its 
hermiticity, the normalization of kinetic terms, etc. We 
again refer to the FeynRules manual for more infor- 
mation [1]. 

After having imported the FeynRules package into 
the current Mathematica session, the model is loaded 
by issuing 



LoadModel[ "filel.fr", "file2.fr", 



] 



3 In this work, we denote as unphysical any field that is not a 
mass eigenstate of the theory. 



where its implementation can be possibly spread among 
the files "filel.fr", "file2.fr", etc., according to 
the convenience of the user. It can also be directly im- 
plemented within the Mathematica session so that 
the function LoadModel is called, in that case, without 
any argument. 

The Feynman rules can be subsequently extracted 
by means of the command 

FeynmanRules [ Lag ] 

where Lag is the Mathematica symbol containing the 
expression of the Lagrangian, written in four-dimensi- 
onal spacetime and employing four-component spinors 
for fermions. The FeynmanRules method extracts the 
interaction vertices included in the Lagrangian Lag so 
that they can be further employed within Mathemat- 
ica for dedicated studies. 

All the interfaces to Monte Carlo event generation 
tools can be invoked in a similar procedure, 

WriteXXX[ Lag ] 

where the sequence of letters XXX takes one of the values 
CHOutput (CalcHep), FeynArtsOutput (FeynArts), 
SHOutput (Sherpa), UFO or WOQutput (Whizard). 
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In the context of supersymmetric theories, the most 
natural and convenient way to construct a Lagrangian 
consists of employing the superspace formalism. There- 
fore, the FeynRules package includes a module allow- 
ing for superfield declarations and Lagrangian imple- 
mentation in terms of superfields. Dedicated functions 
are then provided to convert superfield expressions into 
a form more suitable with respect to the requirements of 
the interfaces to the Monte Carlo event generators [4]. 



3 Implementing mixings in FeynRules 

In this work, we propose an extension of the Feyn- 
Rules package aiming to simplify the declaration of 
the mixing relations linking the unphysical degrees of 
freedom of the theory to the physical fields. This new 
module allows to automatically fill the Definitions 
attribute of the fields, where relevant, and declare the 
mixing matrices as parameters. In addition, we have 
developed an interface generating a CH — h code dedi- 
cated to the diagonalization of the mass matrices of 
the model (see Section 4) after having implemented in 
FeynRules a function allowing for their analytical ex- 
traction from the Lagrangian. In this way, the values of 
all the mixing parameters are derived numerically and 
can be re- imported into FeynRules. 



3.1 Mixing declarations 

For an efficient declaration of the mixing relations, we 
have extended the FeynRules model file structure by 
adjoining a new class dedicated to particle mixings. 
Consequently, all mixing relations among the states can 
be declared on the same spirit as particles, gauge groups 
and parameters, after having been gathered into a list 
dubbed M$MixingsDescription 

M$MixingsDescription = { 

Mix["ll"] == { optionsl }, 
Mix ["12"] == { options2 }, 

} 

Each element of this list consists of an equality ded- 
icated to one specific mixing relation. It associates a 
label, given as a string, ("11", "12", etc.) with a set of 
Mathematica replacement rules defining the mixing 
properties (optionsl, options2, etc.). 

In order to illustrate the choice of options offered to 
the user, we consider the example of Section 2 where 
we have focused on the mixing of the SU(2)l gauge 
bosons. We start by implementing the mixing of the 



W\ and W2 gauge fields, 

, wl - iwl wl + iwt 

W+= " " and W~= " 11 , (1) 

which stems from the diagonalization of the third gen- 
erator of SU(2)l in the adjoint representation. As Eq. 
(1) is purely numerical, i.e., it does not involve any 
model parameter, it can be declared in the model file 
in a very compact form, 

Mix["Wmix"] == { 

MassBasis -> {W, Wbar}, 
GaugeBasis -> {Wi [1] , Wi[2]}, 
Value -> { -[1/Sqrt[2], -I/Sqrt[2]>, 

■Ci/Sqrt [2] , I/Sqrt[2]> } 

} 

The command above declares a mixing relation, dubbed 
Wmix, that can be schematically written as 

MassBasis = Value . GaugeBasis 

where the dot product stands for the usual matrix prod- 
uct. The information on the gauge basis is provided as 
the value of the attribute GaugeBasis which refers here 
to the unphysical fields Wi (Wi [1] ) and W 2 (Wi[2]). 
Similarly, the MassBasis attribute refers to the mass 
basis, containing here the symbols associated with the 
W + (W) and W~ (Wbar) bosons. Finally, the mixing ma- 
trix is given under a numerical form as the argument of 
the attribute Value. 

Some remarks are in order. First, the gauge basis 
only contains unphysical fields, while the mass basis can 
contain either physical fields, unphysical fields or both. 
Particle mixings can therefore be possibly implemented 
in several steps, as illustrated in Section 5.2. Next, spin 
and Lorentz indices can be omitted and the index or- 
dering is defined when declaring the fields (through the 
attribute Indices of the particle class [1]). Finally, if 
some indices are irrelevant, i.e., if they are identical 
for all the involved fields, underscores can be employed 
to simplify the mixing declaration. For instance, the 
three left-handed down-type squarks d c L (sdL [1 , c] ), 
s c L (sdL[2,c]), and b c L (sdL[3,c]) are related to the 
squark gauge-eigenstates Q % [ c (QLs [i , f , c] ), the index 
i being a fundamental SU (2) l index, the index / a fla- 
vor index and the index c a fundamental color index. 
The corresponding declaration reads 

Mix["sdleft"] == { 
MassBasis -> 

{sdL [1 , _] , sdL [2, _] , sdL [3, _] } , 
GaugeBasis -> 

{QLs[2,l,_], QLs[2,2,_], QLs[2,3,_]}, 

} 
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where the mixing matrix is the identity. The underscore 
reflects that the same color index is carried by all fields. 

We now get back to weak gauge boson mixings and 
turn to the neutral sector. We hence focus on the rota- 
tion of the third weak boson W3 and the hypercharge 
gauge boson B to the photon and Z-boson states, 



Wf, 



(2) 



after introducing the a priori unknown weak mixing 
matrix U W A . The computation of the numerical values of 
its matrix elements is addressed by means of the CH — h 
package generated by FeynRules (see Section 4) and 
is only possible if the mixing is declared according to 
the syntax 

Mix["AZmix"] == { 

MassBasis -> {A, Z}, 
GaugeBasis -> {B, Wi[3]>, 
MixingMatrix -> UW, 
BlockName -> WEAKMIX 

} 

The declaration of the gauge and mass bases is simi- 
lar to the case of the charged W bosons, while the at- 
tribute Value has been removed as the numerical value 
of the mixing matrix is not known. The user provides in- 
stead the symbol referring to the mixing matrix (UW) by 
means of the MixingMatrix attribute, without declar- 
ing it as one of the model parameters. This last task is 
internally handled by FeynRules which assumes that 
the mixing matrix is complex and which creates two 
external tensorial parameters, one for the real part and 
one for the imaginary part of the matrix, together with 
one internal tensorial parameter being the matrix itself. 

When a symbol for a mixing matrix is provided, it 
is mandatory to specify, in addition, the name of a Les 
Houches block which will contain the numerical values 
associated with the elements of the matrix. We indeed 
recall that both FeynRules and most of the interfaced 
Monte Carlo event generators order the model param- 
eters according to a structure inspired by the Super- 
symmetry Les Houches Accord (SLHA) [30,31]. In our 
example, we impose the real part of the elements of 
U w to be stored in a Les Houches block WEAKMIX and 
their imaginary part in an automatically created block 
IMWEAKMIX, i.e., a block of the same name with the 
prefix IM appended. 

Implementing model Lagrangians might require to 
explicitly use one or several of the mixing matrices for 



4 Following more standard conventions, the relation of Eq. (2) is 
usually written in terms of the cosine and sine of the electroweak 
mixing angle, as in Section 2. However, we have adopted the 
choice of staying fully general for the sake of the example. 



some of the model interactions, as for the Minimal Su- 
persymmetric Standard Model where the CKM matrix 
is employed in the superpotential [4]. In this case, the 
matrices must be declared according to the standard 
syntax presented in the FeynRules manual, numeri- 
cal values being provided as inputs. This subsequently 
renders the attribute BlockName of the mixing class ob- 
solete and ignored by FeynRules. Contrary, mixing 
matrices automatically declared through a mixing dec- 
laration cannot be employed in Lagrangians. 

In the Standard Model, the CKM matrix Vckm re- 
lates the left-handed down quark gauge-eigenstates dP L 
to the mass-eigenstates 4 as 



d° L — Vckm ■ d-L ■ 



(3) 



To be compliant with the syntax presented so far, a 
symbol for the hermitian-conjugate matrix has to be 
created. To avoid such a complication, the optional at- 
tribute Inverse can be used and set to True, which 
enforces a relation among the mass and gauge bases 
given by 

GaugeBasis = MixingMatrix . MassBasis 



3.2 More advanced cases 

3.2.1 Scalar /pseudo scalar splittings 

When neutral scalar fields are mixing, the gauge eigen- 
states in general split into their real degrees of freedom 
so that one scalar and one pseudoscalar mass basis are 
required. Consequently, a list of two bases is provided as 
argument of the MassBasis attribute, instead of a sin- 
gle basis as in Section 3.1. Consistently, the arguments 
of the attributes Value, BlockName, MixingMatrix and 
Inverse are also upgraded to lists. The first element of 
those lists always refers to the scalar fields, while the 
second one is related to the pseudoscalar fields. It may 
appear that some of the elements of those lists are ir- 
relevant, as for instance when the scalar mixing matrix 
is unknown (MixingMatrix and BlockName are used) 
and the pseudoscalar mixing matrix is known (Value is 
used). The irrelevant list components are in this case 
replaced by underscores, as illustrated with 

Mix["scalar"] == { 

MassBasis -> { {hi, h2}, {al, a2> }, 

GaugeBasis -> { phil, phi2 }, 

BlockName -> { SMIX, _ }, 

MixingMatrix -> { US, _ >, 

Value -> { _, ... } 
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where the (pseudo) scalar mass-eigenstates are repre- 
sented by the symbols hi and h.2 (al and a2). In this 
example, the mixing matrix related to the scalar sec- 
tor is denoted by US and is associated with the Les 
Houches block SMIX. Concerning the pseudoscalar sec- 
tor, a numerical mixing matrix is instead provided (in 
the ellipses). A concrete example is given for the Two- 
Higgs-Doublet Model in Section 5.1. 

3.2.2 Dirac and Weyl fermion mixings 

Several options are left to the user concerning the imple- 
mentation of Dirac fermions mixings. Either one single 
gauge basis is employed, so that FeynRules internally 
takes care of the chirality projectors that appear in the 
related mass terms, or different particle classes can be 
used for the left-handed and right-handed components 
of the fermions. In this case, the GaugeBasis attribute 
refers to a list of two gauge bases instead of to a single 
basis. For both options, the arguments of the attributes 
Value, BlockName, MixingMatrix and Inverse consist 
of lists, the first component being related to the mix- 
ing of the left-handed fermions and the second one to 
the mixing of the right-handed fermions. As for neutral 
scalar mixing, underscores are used for irrelevant list 
elements. 

Lagrangian mass terms for charged Weyl fermions 
are generically written as 




where M stands for the mass matrix and ipi and Xi 
are Weyl fermions which have been assigned an electric 
charge of ±1 for the sake of the example. The diagonal- 
ization of the matrix M proceeds through two unitary 
rotations U and V, 

(4>i\ M"\ fxt\ (xt\ 

: = U : and : = V : , (5) 

W W \xt) \xt) 

which introduces two mass bases. Therefore, all the 
attributes MassBasis, GaugeBasis, Value, MixingMa- 
trix, and BlockName now take lists as arguments (with 
underscores included where relevant). The only extra 
rule to obey to is that the first components of these 
lists are associated with one of the two rotations and 
the second components with the second of them. An 
example is provided in Section 5.3. 



3.3 Vacuum expectation value declarations 

In realistic new physics models, the ground state of the 
theory is non-trivial and fields must be shifted by their 
vacuum expectation value. Since Lorentz invariance and 
electric charge conservation impose that only electri- 
cally neutral scalar fields can get non-vanishing vac- 
uum expectation values, only shifts of (electrically) neu- 
tral scalar fields are allowed to be included in the mix- 
ing relations. This information is encompassed within 
the variable M$vevs which consists of a list of two- 
component elements. The first one refers to an unphys- 
ical field while the second one is the associated vacuum 
expectation value. For instance, the declaration of a 
configuration where two fields phil and phi 2 get non- 
vanishing vacuum expectation values vevl and vev2 
could be performed as 

M$vevs = { { phil, vevl }, { phi2, vev2 } } 

The vacuum expectation values vevl and vev2 must be 
declared as any other model parameter, as described in 
the FeynRules manual [1]. 

3.4 User functions 

Once both the mixing relations and the vacuum ex- 
pectation values have been properly declared, the mass 
matrices of the model can be extracted by means of the 
function ComputeMassMatrix, 

ComputeMassMatrix [ Lag, options ] 

where Lag is the model Lagrangian and the symbol 
options stands for optional arguments. If no option 
is provided, the function calculates all the mass matri- 
ces of the model for which the numerical value of the 
mixing matrix is unknown. It is possible to focus on a 
specific mixing relation whose label is denoted by "11" 
by issuing, in Mathematica, 

ComputeMassMatrix [ Lag, Mix->"11" ] 

For the computation of multiple matrices, the label 
"11" has to be replaced by a list of labels. During the 
computation of the mass matrices, a lot of information 
is by default printed to the screen. This can be avoided 
by including the optional argument ScreenOutput -> 
False in the two command lines above. 

The input information and the result of the Com- 
puteMassMatrix function can be retrieved through the 
intuitive printing functions, MassMatrix, GaugeBasis, 
MassBasis, MixMatrix, BlockName and MatrixSymbol 
which all take as argument the label of a mixing rela- 
tion. A wrapper is also available, 

MixingSummary [ "11" ] 
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which sequentially calls all the printing functions for 
a mixing relation represented by the label "11" and 
organizes the output in a readable form. 

The FeynRules method to extract analytically a 
mass matrix is fully generic and can be employed to 
compute any matrix M defined by the Lagrangian 

C masa = B\ M B 1 , (6) 

where B\ and B2 stand for two field bases possibly dif- 
ferent. The calculation of the matrix M is achieved by 
issuing 

ComputeMassMatrix [ Lag, 

Basisl -> bl, Basis2 -> b2 ] 

where the symbols bl and b2 are associated with the 
bases B\ and B2 and refer to lists of fields. In this case, 
the printing functions introduced above are not avail- 
able. 

4 Automated spectrum generation 

4.1 The ASperGe package 

The computation of the unknown mixing matrices nec- 
essary for diagonalizing all the model mass matrices 
can in general only be achieved numerically. To this 
end, we have developed the CH — h program ASperGe. 
It includes a set of CH — h source files (stored in the 
subdirectory src), coming together with the related 
header files (stored in the subdirectory inc), that can 
be split into model-independent and model-dependent 
files. For an efficient use of the ASperGe program, 
it has been entirely embedded within the FeynRules 
package. Therefore, only a brief discussion of the struc- 
ture of the code is presented in this paper. More in- 
formation, such as a doxygen documentation, can be 
found on the ASperGe webpage [32]. 

The set of model-independent files contains, on the 
one hand, several tools dedicated to matrices and their 
diagonalization (MassMatrix . cpp, MassMatrix.hpp as 
well as Matrix . hpp) . On the other hand, the ASperGe 
code is based on an internal format for parameters, de- 
fined in the source files Par. cpp, CPar.cpp, RPar.cpp 
and in the associated header files. This format is in- 
spired from a SLHA structure and the corresponding 
mapping is encoded into the files ParSLHA . cpp, SLHA- 
Block.cpp, and in the associated header files. Finally, 
printing and string manipulation routines are included 
in the files tools . cpp and tools . hpp and the program 
conies with a makefile. 

All the model dependency is included in the two files 
Parameters . cpp and Parameters . hpp as well as in the 
core program implemented in the main, cpp file. 



The information encompassed in the two parameter 
files is threefold. First, the SLHA structure ordering 
the external parameters is encoded in terms of blocks 
and counters. Next, the definitions of the internal pa- 
rameters as functions of the other model parameters 
are implemented, where a proper running of the AS- 
perGe program is only guaranteed if the parameters 
do not depend on the masses and mixing matrices to be 
computed. Finally, the analytical formulas of the mass 
matrices to diagonalize are included. 

The main program (main . cpp) starts with the dec- 
laration of the different mass matrices of the model. 
Links to the relevant elements of the mass basis are then 
implemented by means of the associated PDG codes, 
which allow to assign the mass eigenvalues to each of 
the physical particles, the ordering of the PDG codes 
following the mass ordering. 

4.2 Interfacing the ASperGe package to FeynRules 

The ASperGe package can be entirely generated, for 
a given particle physics model, from the FeynRules 
model information by means of a dedicated interface 
which works as for the other FeynRules interfaces. It 
is then called by typing, in a Mathematica session, 

WriteASperGe [ Lag, Output -> dirname ] 

where the symbol Lag stands for the model Lagrangian 
and Output ->dirname for an optional argument indi- 
cating the name of the directory where to store all the 
created files. If unspecified, the directory ModelName_MD 
is employed, ModelName being the name of the Feyn- 
Rules model. 

The interface first extracts all the relevant mass ma- 
trices from the Lagrangian Lag by means of the func- 
tion ComputeMassMatrix introduced in Section 3.4. It 
then writes, in addition to model-independent files de- 
scribed in Section 4.1, the three model-dependent files 
main. cpp, Parameters . cpp and Parameters . hpp, to- 
gether with one data file Externals . dat (stored in the 
subdirectory input). This last file contains the numeri- 
cal values of the external parameters of the model, nec- 
essary for the numerical evaluation of the mass matri- 
ces. When running the code (see Section 4.3), the user 
can update this file or even employ a different file ac- 
cording to his needs. 

The numerical matrix diagonalization performed by 
ASperGe is based on Gsl functions relying on the her- 
miticity of the mass matrices which employs symmetric 
bi-diagonalization followed by QR reduction. This con- 
trasts with existing diagonalization packages developed 
in the framework of FeynArts [33] and CalcHep [34] 
that are based on Jacobi-type iterative algorithms. A 
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hermiticity check is therefore performed by the inter- 
face before writing down the output. Since the mass 
matrix Af related to charged fermions is by construc- 
tion non-hermitian, the matrices Aft Af and Af Aft are 
employed instead, which allows to obtain left-handed 
and right-handed fermion mixing matrices separately. 

It is also possible to focus on one or several specific 
mixing relations. In this case, the Mix option, already 
introduced in the context of the ComputeMassMatrix 
function, has to be used, 

WriteASperGe [ Lag, Mix -> {"11", "12"} ] 
We refer to Section 3.4 for more information. 

4.3 Running ASperGe 

Since the ASperGe package is based on Gsl functions, 
it is mandatory to have the Gsl libraries installed on 
the system. Then, if the g++ compiler is available, the 
makefile generated by FeynRules can be employed di- 
rectly. Otherwise, it must be first edited accordingly to 
include proper compiler information. 

Once compiled, ASperGe can be executed by typ- 
ing in a shell 

./ASperGe <infile> <outfile> 

where the arguments indicate in which file the numeri- 
cal value of the external parameters must be read (<in- 
f ile>) and where to store the output file (<outf ile>). 
This file contains, in addition to the input parameters, 
the computed numerical values of the mixing matrices, 
split in terms of their real and imaginary parts accord- 
ing to the SLHA conventions, as well as all the masses 
of the physical states (stored in the SLHA block MASS). 
In order to execute ASperGe with all the default set- 
tings as generated by the FeynRules interface, it is 
sufficient to type in a shell 

./ASperGe input/externals . dat output/out . dat 

Both the compilation and the execution of the pro- 
gram can be performed from the Mathematica ses- 
sion, by issuing 

RunASperGe[ ] 

This also loads the SLHA parameter file out . dat back 
into the FeynRules session, so that it can be fur- 
ther employed, e.g., to generate a UFO model. Infor- 
mation about the run of ASperGe can be found in 
the file ASperGe . log stored in the same folder as the 
executable. 

It is also possible to diagonalize specific mass ma- 
trices of the model by executing 

./ASperGe <infile> <outfile> ml m2 ... 

where ml, m2, etc., are the names of the mixing matrices 
to be computed. 



5 Illustrative examples 

In this section, we illustrate the features of the AS- 
perGe program and its interface to FeynRules by 
choosing three extensions of the Standard Model with 
non-trivial mixing relations, i.e., the Two-Higgs-Dou- 
blet Model (2HDM), the Minimal Left-Right Symmet- 
ric Standard Model (LRSM) and the Minimal Super- 
symmetric Standard Model (MSSM). We modify their 
original FeynRules implementations 5 [2,4] to accom- 
modate for the mixings as described in Section 3. We 
then employ the ASperGe program (see Section 4) to 
numerically calculate some of the mass and mixing ma- 
trices of these models. 

5.1 The general Two-Higgs-Doublet Model 

The 2HDM is one of the simplest extensions of the 
Standard Model, with respect to which it only con- 
tains a second weak doublet of scalar fields. Following 
the conventions of the original FeynRules implemen- 
tation [2], both Higgs fields fa and fa carry the same 
hypercharge so that they can always be redefined by 
means of U(2) transformations [35,36,37,38]. Adopting 
the so-called Higgs-basis, the two doublets read 

where only the neutral component of the fa field ac- 
quires a vacuum expectation value v. Moreover, the 
Goldstone bosons G ± and G° as well as the charged 
Higgs field 77 ± are not required to be further rotated, 
so that only the mass matrix of the neutral fields H°, 
R° and 1° must still be diagonalized. 

We extend the 2HDM FeynRules implementation 
described in Ref. [2] by first indicating that the sec- 
ond component of the fa field, represented by the sym- 
bol phil, acquires a no n- vanishing vacuum expectation 
value v labeled by the symbol vev, 

M$vevs = { {phil [2], vev} } 

as shown in Section 3.3. Then, we choose to implement 
the mixing of the H°, R° and 1° fields to the physical 
hi, hi and /13 fields in a two-step manner. In a first 
stage, the gauge eigenstates are split into their scalar 
and pseudoscalar components, 

Mix["lp"] == { 

MassBasis -> { {HO}, {GO} }, 
GaugeBasis -> { phil [2] }, 

5 Since no previous implementation of the LRSM exists, we take 
the opportunity to provide the relevant details in Section 5.2. 
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Value 



-> { {{1» > 



Mix["2p"] 
MassBasis 
GaugeBasis 
Value 



{ 



-> { -CRO}, {10} }, 

-> { phi2[2] }, 

-> { ill}} } 



following the syntax introduced in Section 3.1 and Sec- 
tion 3.2 and making use of the self-explained symbols HO 
and R0 (GO and 10) for representing the (pseudo) scalar 
degrees of freedom. Similarly, we can employ the mix- 
ing infrastructure to map the charged components of <pi 
and 4>2 to the physical fields G + and H + by means of 
a 1 x 1 identity matrix. Since this procedure is trivial, 
we omit any further details from the present manuscript 
and refer to the model implementation [32] . In a second 
stage, the rotation to the physical fields, represented by 
the symbols hi, h.2 and h.3, is declared as 

Mix["ls"] == { 

MassBasis -> { hi, h2, h3 >, 
GaugeBasis -> { HO, R0, 10 }, 
MixingMatrix -> NH, 
BlockName -> NHMIX 

} 

where we associate the symbol NH to the correspond- 
ing mixing matrix and assign the Les Houches block 
(IM)NHMIX to the numerical value of its elements. The 
neutral squared mass matrix M 2 can then be derived 
from the model Lagrangian (represented by the symbol 
L2HDM) by typing, in the Mathematica session, 

ComputeMassMatrix[L2HDM, Mix->"ls"] 

As a result, one recovers the well-known expression de- 
pending on the most general scalar potential param- 
eters Xi and fii (see Ref. [2] for further information), 



M 2 



( 2X lV 2 



%t[X 6 ]v 2 m 

,,2 



■+A E 



\-%{\ 6 }v 2 





— -A 5 



(8) 



v 2 - J 



after having introduced the squared mass of the charged 
Higgs boson m| = l/2A3W 2 +^2 and removed two of the 
/i-parameters by means of the potential minimization 
conditions, [i\ = — Aiv 2 and ^ = —l/2Xev 2 . 

The numerical value of the unitary matrix U diago- 
nalizing M 2 is obtained by generating and making use 
of the ASperGe package, as shown in Section 4. We 
fix, adopting a representative benchmark scenario, the 
Higgs potential parameters to Ai = A2 = A3 = 1.0, A4 = 



Table 1 Field content of the LRSM, given together with their 
representation under the SU(3) C X SU(2) L X 5(7(2) R X U(1) B - L 
gauge group. The SU(2) L = 1, 2) and SU(2) R = 1, 2) 

fundamental index structure is explicitly indicated. 



Field 



Components 



Representation 



Q'l 

Qui' 

1 1 
L L 

Lav 



(3,2,1,|) 
(3,1,2*,— ; 
(1,2,1,-1) 
(1,1,2*, 1) 



1 l' 



<p0 <p+ 

<p''- <p'0 



*1 + 



-k=A + A ++ 

AO 1 A + 

^R V2 R 



(1,2,2*,0) 

(1,3,1,2) 

(1,1,3,2) 



0.5, A 5 = 0.4, A 6 = 0.3, A 7 = 0.2 and ^ = 6 • 10 4 GeV. 
This leads to the three mass eigenvalues 



(m hl ,m h2 ,m h3 ) = (285,327,379) GeV, 

whereas the mixing matrix reads 

-i N 

U= I 0.784 -0.621 
-0.621 -0.784 0, 



(9) 



(10) 



These results are in good agreement with those ob- 
tained by means of the TwoHiggsCalc calculator [16]. 



5.2 The Minimal Left-Right Symmetric Standard 
Model 

The LRSM [39,40,41,42,43,44,45] is an extension of 
the Standard Model with an enlarged SU(3) c xSU(2) L x 
SU(2)r x U(1)b-l gauge symmetry. In this model, the 
fermionic degrees of freedom of the Standard Model 
lying in the trivial representation of SU(2)l are col- 
lected into SU(2)r doublets, as shown in the first part 
of Table 1 where the model matter field content is pre- 
sented together with the associated quantum numbers. 
In addition, the symmetry-breaking mechanism down 
to electromagnetism is also more involved, relying on 
an enriched Higgs sector (see the second part of the 
table). 

The LRSM Lagrangian consists of standard kinetic 
and gauge interaction terms for all fields as well as of 
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the Yukawa interactions 

+ QrY^^Ql + L^&Ll (11) 
+ L c L yf } A L L L + L R yf ] A R L R + h.c. . 

In this equation, all indices are understood, the matri- 
ces yq and are 3x3 matrices in flavor space and 
the superscript c indicates charge conjugation 6 . More- 
over, gauge invariance is ensured by the introduction of 



the hatted fields 

£/ = tijj'i'&r , L Ll = ajLi , L R = e i,j 'L Rf , (12) 

where the rank-two antisymmetric tensors with lower 
and upper indices are defined by e±2 = — e 12 = 1. Intro- 
ducing Higgs mass parameters \ii and quartic interac- 
tion strengths A,, pi and a.i, the scalar potential reads 



£ H = M 2 Tr[<?t<?] - Ai (Tr[&$]) - \ 2 Tr[&$&$\ - -A 3 (*Tr[^*] + Tr^t - A 4 Tr[<P t ^ t <? t *] 

- -A 5 (Tr[#<2>*] - Tr[# t # t *]) - -X 6 (Tr[^S^S^] + Tr[<P*<2><£*<2>]) + p\ {jx\A\A l \ + Tr[A R A R ]j 

- Pl {tt[AIA l } 2 + Tr[A R A R ] 2 ) - P2 (tt[A[A l AIA l }+Tt[A r A r A^A r }) - p 3 Tr[A{A L ]Tr[A R A R ] (13) 

- ai Tr[&$](Tr[AlA L } + Tr[^^ R ]) - a 2 (Tr[A R &$A R ] + Tr[A\$& A L ]} 

- ^(^[A^^Ar] + TtIAIS^Al]) . 

I 



Since the corresponding FeynRules model description 
is standard, we refer to the FeynRules manual [1] and 
leave all implementation details out of this work. 

In the LRSM, the symmetry-breaking mechanism is 
performed in two steps. At high energy, the SU(2)l x 
SU(2) R x U(1)b-l gauge symmetry is spontaneously 
broken to the electroweak symmetry, the latter being 
subsequently broken to electromagnetism at a lower 
scale. Consequently, the neutral components of the sca- 
lar fields get vacuum expectation values at the mini- 
mum of the potential, (<£°) = v/y/2, (<P'°) = v'/-/2 and 
(A° L R ) = vl,r/V2, by which they are shifted. In the 
rest of this section, we focus on the mixing of the neu- 
tral Higgs fields and illustrate the way to implement a 
two-stage field rotation. For all the other mixing rela- 
tions of the LRSM, we refer to the implementation [32]. 

We first assume, motivated by neutrino mass and 
kaon system data [45,46], that vl = v' « 0. Next, we 
implement the rotation associated with the diagonal- 
ization of the third generator of SU(2) in the adjoint 
representation as 

Mix["2a"] == { 

MassBasis -> { DLpp, DLO }, 

GaugeBasis -> { DL[1], DL [2] >, 

Value -> { {1/Sqrt [2] , -I/Sqrt[2]}, 



6 We recall that the components of the field Qr are charge con- 
jugate (see Table 1). 



-Cl/Sqrt[2] , I/Sqrt[2]» 

} 

depicting the example of the SU (2)l Higgs triplet. The- 
se replacement rules translate the rotation of the A\ 
and A\ fields, represented by the DL[1] and DL[2] 
symbols, to the A and A ++ states labeled by DLO and 
DLpp. Then, the neutral fields A° L , A° R , <P° and <P'° , 
represented by the symbols DLO, DRO, phi [1,1] and 
phi [2,2], mix to four scalar degrees of freedom h\, h%, 
and h®, two physical pseudoscalar Higgs bosons a? 
and a° and two Goldstone bosons G® and G° to be 
eaten by the Z and Z' vector fields when getting mas- 
sive. Introducing the corresponding symbols hOl, h02, 
h03, h04, aOl, a02, G01 and G02, these rotations are 
implemented as 

Mix["2e"] == i 

MassBasis -> { {hOl ,h02 ,h03 ,h04} , 

{G01,G02,a01,a02> }, 

GaugeBasis -> 

{ DLO, DRO, phi [1,1] ,phi[2,2] }, 
MixingMatrix -> { UHN.UAN }, 
BlockName -> { HMIX,AMIX> 

> 

The two symbols UHN and AUN respectively denote the 
scalar and pseudoscalar mixing matrices, the numerical 
value of their elements being included in the two Les 
Houches blocks (IM)HMIX and (IM)AMIX. 
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Typing, in Mathematica, the commands 

ComputeMassMatrix[Lag, Mix->"2e"] 
MassMatrix["2e" , "S"] 

allows to calculate both the scalar and pseudoscalar 
mass matrices and display the scalar squared mass ma- 
trix Ai 2 to the screen. It reads, 



( A 








°\ 





B (a i - 




v v R 


(ai - 


\- a 3 )v v R 


C 





\o 








DJ 



Ai' 



where we have introduced the quantities 
A = ^(ai + a 3 )v 2 - p\ + ^p 3 v 2 R , 
B= ^(ai + a 3 )v 2 - fij + 3(pi - p 2 )v R , 



(14) 



C = 3(Ai + \ 2 )v 2 + + a 3 )v 2 R - p\ , 



(15) 



D = (Ai+4A 3 + A 4 + A 6 )w 2 + -{a 1 +a 2 ) 



a 2 )v 2 R 



Mi 



In order to numerically compute the unitary matrix 
U diagonalizing Ai 2 , we use the ASperGe program, 
generated from FeynRules by issuing 

WriteAsperge [Lag] 

after fixing the external Lagrangian parameters of the 
model to Ai = A2 = A3 = A4 = A 6 = 0.1, ot\ = 0.1, 
a 2 = 0.3, a 3 = 0.1, p\ = 0.1, p 2 = and p 3 = 0.5. 
In addition, the //-terms are deduced from the mini- 
mization conditions of the scalar potential and the vac- 
uum expectation values are taken as v = 248 GeV and 
v R = 6000 GeV. It should be noted that the relevance 
of such numerical values is going beyond the scope of 
this paper, and could be addressed by means of external 
packages such as the one presented in Ref. [47]. Once 
ASperGe is executed, one obtains 



U = 



/0 -0.041 0.99 0\ 
1 
10 

\0 0.99 0.041 0/ 



(16) 



the mass eigenvalues being 

(mi, m 2l m 3 , m 4 ) = (111, 1905, 2324, 2686) GeV. (17) 

5.3 The Minimal Supersymmetric Standard Model 

In supersymmetric extensions of the Standard Model, 
each of the model's degrees of freedom comes accom- 
panied by a superpartner with opposite statistics. The 
minimal version of such theories, the so-called MSSM 



[48,49], has been originally implemented in FeynRules 
by making use of its superspace module [4]. We hence 
refer, on the one hand, to Ref. [4] for notations and 
conventions and, on the other hand, to the new model 
implementation [32] for more information on the way 
in which particle mixings have been implemented. In 
the rest of this subsection, we employ the MSSM im- 
plementation to illustrate fermion mixing declaration. 

After electroweak symmetry breaking, the left-hand- 
ed, two-component Weyl fermionic, gaugino and hig- 
gsino fields W and HJ / mix to the chargino eigen- 
states x ± ■ Introducing the two rotation matrices U and 
V (labeled by the symbols UU and VV), this mixing is 
declared through an instance of the mixing class, 



Mix["3d"]=={ 
MassBasis 

GaugeBasis 

BlockName 



-> { {chmw[l] ,chmw[2]>, 
{chpw[l] ,chpw[2]} }, 

-> { {wowm , hdw [2] } , 
{wowp,huw[l] } }, 

-> -CUMIX.VMIX}, 



MixingMatrix -> {UU.VV} 

} 

In this list of replacement rules, wowm (hdw) and wowp 
(huw) are the labels of the negatively and positively 
charged wino (higgsino) states, the related mass eigen- 
states being represented by the symbols chmw and chpw. 
In addition, we link the mixing matrices U and V to the 
Les Houches blocks UMIX and VMIX. 

Like in the previous subsections, we employ Feyn- 
Rules to extract the corresponding tree-level mass ma- 
trix Ai' from the Lagrangian, 

M 2 \/2mw sin (3^ 
y \f2m w cos /3 p 



M' 



(18) 



where mw denotes the TV-boson mass, p the superpo- 
tential Higgs mixing parameter, M 2 the supersymmetry- 
breaking wino mass and tan (3 is defined as the ratio of 
the two neutral Higgs field vacuum expectation values 
tan/3 = {H°/H°).^ 

As in the original implementation, we choose the 
typical minimal supergravity point SPS la [50] as a 
benchmark scenario. Then, we generate the MSSM AS- 
perGe program and use it to compute numerically the 
mixing matrices U and V, 



U = 



0.918 -0.397 
-0.397 -0.918 



0.974 -0.226\ 
-0.226 -0.974 ) ' { ' 



as well as the corresponding mass eigenvalues 

(m x ±,m x ±) = (176,382) GeV. (20) 

Those values are in good agreement with those returned 
by commonly used MSSM spectrum generators. 
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6 Summary 

In this paper, we have presented an extension of the 
FeynRules package dedicated to the automated gen- 
eration of the particle mass spectrum and mixing struc- 
ture associated to any Lagrangian-based quantum field 
theory. The new module is based on the introduction of 
a new structure for particle mixing declaration allow- 
ing, on the one hand, for the analytical computation of 
all the model mass matrices, and, on the other hand, 
for the generation of a CH — h program dubbed AS- 
perGe, yielding the numerical evaluation of the associ- 
ated rotation matrices. We illustrate the strength of this 
new module in the context of the Two-Higgs-Doublet 
Model, the Minimal Left-Right Symmetric Standard 
Model and the Minimal Supersymmetric Standard Mo- 
del. 
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